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Abstract 



We derive renormalization group equations for neutrino masses, leptonic 
mixing angles and CP-violating phases running at energies above the seesaw 
threshold, both in the Standard Model and in the Minimal Supersymmetric 
Standard Model extended with three right-handed neutrinos. With these 
equations, we carry out a systematic study of the radiative correction that 
may arise to neutrino parameters, via their renormalization group evolution 
from the electroweak scale (Mz) to the scale of Grand Unified theories (Aqut)- 
We study in detail three typically interesting neutrino mass patterns: normal 
hierarchy, near degeneracy and inverted hierarchy. Magnitudes of possible 
corrections in each case are carefully investigated. We also emphasize the 
significance of CP-violating phases in controlling the evolution behavior of all 
neutrino parameters. 
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I. INTRODUCTION 



Experimental information on neutrino masses and mixings has opened up a new playground 
for efforts [1] of a better understanding of the Yukawa couplings in the Standard Model 
(SM) and its extensions. In the lepton sector, it is conceptually natural to use the elegant 
seesaw mechanism [2-4] to give masses to the SM neutrinos. In the most simple version of 
this mechanism, the dimension 5 neutrino mass term is the low energy relic of some more 
fundamental theories with very heavy right-handed neutrinos. So in this framework, we 
often need to relate physics at vastly different energy scales. 

The only way to compare the high energy theoretical prediction and the low energy ex- 
perimental observation is to use renormalization group equations (RGEs). It has been found 
that corrections arising from renormalization group (RG) evolution can be very significant 
for leptonic mixing angles and neutrino mass splitting, especially in the case of nearly degen- 
erate left-handed neutrinos. So in principle, the RG correction should not be neglected in 
the discussion of models suggested at high energy scales. With RGEs derived in Refs. [5-7] , 
early discussions [8-16] of the RG running effect are mainly concerned with neutrino mass 
(or Yukawa coupling) matrices. The evolution of neutrino masses, leptonic mixing angles 
and possible CP-violating phases is studied by diagonalizing the relevant Yukawa matrices 
at different energy scales, and the behaviors are discussed in a numerical or semi- analytical 
way. However, authors of Refs. [17-20] have also emphasized the significance of RGEs for 
individual neutrino parameters. Such equations not only make it possible to predict the 
evolution behavior of each parameter [21-24], but also can help people appreciate inter- 
esting features such as the existence of (pseudo-) fixed points in the evolution of mixing 
angles and CP-violating phases [14,16,17,25,26]. The derivation of these equations below 
the seesaw threshold has been done in Refs. [17,18,27,23]. And based on such equations, a 
comprehensive study of the RG evolution of neutrino parameters from the electroweak scale 
to the seesaw threshold has been carried out in Ref. [27]. 

However, RG corrections above the seesaw threshold sometimes are as important as or 
even more significant than those below the threshold [8,15,28-30]. Since the physics re- 
sponsible for neutrino mass generation is more likely to exist at the scale of Grand Unified 
theories, a systematic study of the RG correction above the seesaw threshold should be nec- 
essary. And this is one of our main concerns in this work. In much the same spirit as of Ref. 
[27], we derive RGEs for individual neutrino parameters running above the seesaw thresh- 
old, under the condition that eigenvalues of the Yukawa coupling matrix that connects left- 
and right-handed neutrinos are hierarchical. The contribution from the largest eigenvalue 
of this matrix is explicitly shown. By setting this contribution to zero, we can regain RGEs 
obtained earlier in Refs. [27,23], which are valid at energies below the seesaw threshold. 

The second purpose of this work is to carry out a systematic study of the radiative 
correction that may arise via the RG evolution in the full energy range from the electroweak 
scale (Mz) to the scale of Grand Unified theories (Agut)- To demonstrate main features 
of possible corrections, we study in detail three typically allowed neutrino mass patterns: 
normal hierarchy, near degeneracy and inverted hierarchy. 

The paper is organized as follows. We write down full one-loop RGEs for individual 
neutrino masses, leptonic mixing angles and CP violating phases in Section II, with a brief 
discussion. Then in Section HI, we carry out a systematic study of the correction that may 
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arise during the RG evolution from Mz to Aqut- in theories with the three neutrino mass 
patterns mentioned above. Section IV is devoted to a Summary. 

II. ANALYTICAL FORMULAE FOR NEUTRINO PARAMETERS RUNNING 

ABOVE THE SEESAW THRESHOLD 

Extended with three right-handed neutrinos, the Lagrangian giving mass to leptons in the 
SM is 

-£i = ElHYil^ + 'WlH'^Y.u^ + ]^14^MrUr + /i.e.; (1) 

and in the MSSM is 

-Ci = E^HiYilji + EEH2Y,iyji + ]^^MrVr + h.c, (2) 

where, E^, and denote 5'C/(2)i-doublets, right-handed charged leptons and right- 
handed neutrinos, respectively. Both in Eqs.(l) and (2), the scale of Mr is expected to be ex- 
tremely high, since there is not a protective symmetry. Around this energy scale, mass eigen- 
states of Mr are successively integrated out at their respective masses (Mi < M2 < M3), 
giving rise to a series of effective theories at different energy scales [28] . Then at energies 
below the lightest right-handed neutrino mass, we obtain the dimension 5 effective mass 
term for left-handed neutrinos 

= -]^L^ ■ n ■ ^^El + h.c, (3) 

where k, is the effective Yukawa coupling matrix, and $ is H'^ in the SM but is H2 in the 
MSSM. Since k, is calculated from Y^, and Mr by decouphng right-handed neutrinos at 

successive energy scales, step by step, the relation between n and Yi,, Mr is complicated. 
Only in the most simplified procedure when all right-handed neutrinos are decoupled at 
a common scale, can we obtain (at that chosen scale) a simple equation because of the 
tree-level matching condition [28]: 

K = Y^Mi'Yj . (4) 

Then, when the Higgs field acquires a non-zero vacuum expectation value ($) = {v 0)^ 
during the electroweak symmetry breaking, Eq.(3) yields an effective mass matrix M^, = v'^k 
for left-handed neutrinos. In the SM, v ^ 174 GeV; and in the MSSM, v ~ 174 sin /5 GeV. 

Since right-handed neutrinos are to be decoupled at their respective thresholds, it will be 
too complicated to derive RGEs for neutrino parameters between these thresholds. So we 
shall be less ambitious than solving the whole problem, but shall simplify it by (a) decoupling 
all right-handed neutrinos at a common scale, which we take to be M3, and (b) limiting our 
derivation only to the case when eigenvalues of Yi, are hierarchical. Here, assumption (a) 
can be justified when the RG evolution through right-handed neutrino thresholds is not too 
dramatic (As has been demonstrated in Refs. [15,28], this is not always the case). And 
assumption (b) is also well motivated since in a large class of high energy models considered 
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in literature (such as those based on U{1) symmetry or those base on S'O(IO) Grand Unified 
Theories), there is often a certain similarity or even identification between and the up- 
quark Yukawa coupling matrix. Such a similarity of matrices should lead to some likeness 
between their eigenvalues. 

One-loop RGEs for k, running from Mz to Mi and those for Y^y and Mr running through 
right-handed neutrino thresholds to Aqut have been given in Ref. [28]. For readers' con- 
venience, we have collected a part of them together with those for Yi in Appendix A. To 
discuss the KG evolution of neutrino parameters in the full energy range from Mz to Aqut, 
it is convenient to make use of k, also at energies above the seesaw threshold [24,30]. In this 
energy range, we find from Eqs.(4), (A. 6) and (A. 7) 



at 



(5) 



where i = In/x with being the energy scale, and details of and Ni^ are given in Eq.(A.15) 
in the Appendix. 

For neutrino masses, leptonic mixing angles and CP-violating phases, k is diagonalized 
by a unitary matrix t/^: 



K = Ui^k'U'^; k' = diag {ki, k2, ka} , 



(6) 



where ki (for i = 1,2, 3) at Mz are proportional to left-handed neutrino masses : ki (Mz) = 
rrii (Mz) jv^. At energy scales below Mi, Yi is always diagonal during the RG evolution, if 
it is diagonal at the beginning. In such a basis, the leptonic mixing matrix is ?7mns = ^k- 
However, Yi can not be kept diagonal above the seesaw threshold when Eq.(A.5) is used 
for its RG evolution. In this case, there is the contribution to t/MNS from diagonalizing 



Hi = Y^Y}: 



Hi = UiH[U}; H', = diag {yl, y^} , 
For the MNS matrix, a convenient parametrization can be found in Ref. [31] 



(7) 
(8) 



MNS 



/ 



-iS 



CyS^C C^SyS^ 



CxCySz 



C^CyC S^SyS^ 



-CxSyC 



CySxSz 



Crr 



C^Sy 
CyCz 



\ 



1/ 



(9) 



where Cx = cos^a;,Sa; = sin^j; and so on. ^The merit of this parametrization is that the 
Dirac phase 5 does not appear in the neutrinoless double beta decay, while Majorana phases 



^Comparing to another parametrization used in the Hterature [32], apart from the notation of 
mixing angles 6*12 6x, 623 <^ &y and 6*13 <^ 9z, the main difference resides only in the Majorana 
phases: ai <^ 2(p — 6), a2 2{a — 6), while the Dirac phase 6 is the same. 
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p and a do not contribute to the leptonic CP violation in neutrino oscillations. Thus two 
different types of phases can be separately studied in different types of experiments. 

At energies below the seesaw threshold, RGEs for the running of left-handed neutrino 
masses (mi, 7712,^3), leptonic mixing angles {0x,9y,9z) and CP-violating phases {6,p,a) 
have been derived in Refs. [17,18,27,23]. In order to obtain the same kind of formulae at 
energies above the seesaw threshold, we also need a parametrization of Y,^. We find that 
the derivation is most straightforward if Yi, is parameterized in the diagonal basis of k by 
(h^^Y.YJ): 



= U^KUI- K^yl- diag {r?, rj, l} , 



(10) 



where (with ci = cos^i, Si = sin^i and so on) 



/e'"^! \ /I 



v 



Ci Si 
V -Si Ci / 



C2 



\-S2 




(11) 



As mentioned above, we shall concentrate on cases in which eigenvahies ofY,, are hierarchical, 
i.e. << r| << 1. So the contribution of rf and r| to the running of neutrino parameters 



can always be neglected. This is equivalent to taking rf 



when deriving RGEs. 



From Eqs.(lO) and (11), it is obvious that only 6'i, 6*2,^1 and 02 in U^, contribute to the 
evolution of parameters mi, m2, ms, 9x,9y,0z,S, p and a. 

Before writing down all the analytical formulae, we remark that RG corrections to mixing 
angles and CP- violating phases come from three separable sources. To clarify this point, 
we need Eq.(A.15): A^^ = C^Hi + C^H^. In Eq.(5), only A^^ contributes to the evolution of 
mixing angles and CP-violating phases. So each of the two terms in A^^ is a source of RG 
corrections. Also, there is a contribution from diagonalizing Yi, and this is the third one. 



The contribution from C]^Hi in A^^ is proportional to C^y^ (we have omitted the 
contributions from y1 and for obvious reasons). This contribution is exactly the 
same as that governs the evolution of neutrino parameters at energies below the seesaw 
threshold. Analytical formulae of this contribution have been derived and extensively 
discussed in Refs. [27,23]. 

The contribution from C^i^^ in A^^ is proportional to C^yl- While y-r can be of 0{1) 
only in the MSSM when tan /3 is large, it is quite natural for yi, to be of the same 
magnitude as the top Yukawa coupling. Furthermore, just like the contribution from 
C^Hi, the contribution from C^H^ can also be resonantly enhanced when eigenvalues 
of K are nearly degenerate. This is because both contributions contain such enhancing 
factors as C,~j^ (for i < j; i,j = 1, 2, 3). Note that 



ki kj 
ki -\- kj 



= 1,2,3. 



(12) 



• The third contribution comes from diagonalizing Yi and is proportional to C^y^. Dif- 
ferent from that of C^H^, there are no enhancing factors in this contribution other 
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than functions of mixing angles, such as ^ = (sin ^ or s^. ^ = (sin^^;) ^ etc., which 
appear mostly in RGEs of CP- violating phases and are important only when (or 
etc.) « 0(1). 

Now following the same procedure as described in Refs. [18,27] but including contribu- 
tions from and Yi (yet with the above explained simplifications), we obtain the following 
full one-loop RGEs for left-handed neutrino masses {1711,1712,171^), leptonic mixing angles 
{9x, Oy, 6z) and CP-violating phases {5, p, a) running at energies above the seesaw threshold. 

^For the running of left-handed neutrino masses (at : TTij = V ki ] i — 1, 2, 3): 

ki — (liki , 2 = 1, 2, 3 ; 

ai = a^ + 2Clyl (s^sj - 2csCa,CyS^SySz + clclsfj + 2C'^ylsl ; 

, Q W 2/22,0 _l22 2\, o/^iv 2 2 2 

02 = 0;^ -h '^'^^Vt [Cr^Sy -\- zcgC^CyS^SySz + CyS^s^ J + zo^y^CgSj^ ; 

03 = a« + ^Cly^ycl + 2C'',ylc\cl . (13) 

For the running of leptonic mixing angles and CP- violating phases (c^ = cos 5, = 
sin 5, C(^_p) = cos {5 — p) , S(5_p) = sin {5 — p) and so on): 



0. = C\yl \ ^ [c(,-.)C... [si - cy:) - ( 



12 



Cl3 



{p{8—p)^xSy CpCxCySz^ Cl3 ■ Sp {s(^g_p^SxSy + SpCxCyS- 



,23 



'^{&-(j)'^xSy + C^CyS^Sz 



C23 ■ ( 



CySxSz 



CxCySz 



+c:yi 



C12 

7 C,i3 S S^^ 

S13 



C12 • S(p_.)S(^,_^2) C2S1S2 



CiC2S2SxSz , I C^C^2 \ CiC^CxSiSz 



+ 



+C\yl 



C\C2Sz 



-C(p-a+<Ai-<A2)'^2-SlS2 (4 - 4) } ; 



- C23 • s^s 

/ 2 2 _ 2\ 
^2 ) 



(14) 



Oy = 



+ 



Cl3 



,23 



{c{S-p)^xSy — CpCxCySz^ + Cl3 • •5(5-p) {^{5-p)^xSy + SpCxCyS^ 
{'^(S-a)^xSy + C^CySxSz^ + (23 ' ^(S-<j) {^{S-(7)(^xSy — S^CySxS^ Ca;Cy| 



dF 

^Throughout this work, for what ever F : F = 167r^— t = Ina, with u being the energy scale. 

at 



6 



+ Cl3-S(,_,)S^ 1 — ^ 



+ <,23 ■ 



L V U3 I \ (,23 

+2c5Ca;C^ {cls\ - S2) + 2 (c(^+^^)Ca;S2 + C^^^+^^)C2SlS^ CiC^CyCzSyS 



+ (c(5-p-0i)S2Sa; - Ci^S-<T-<t>2)^'2CxSl 
+2 -C(p_CT+0^_^2)'^3;Sa;Sj/ (l + 



C1C2 



+ (c, 



(5+p-o-+<;/)i -^2) ^a: C(<5-p+(7-(/)i +</>2) *a: , 



) CyS^ C2CySiS2^ ; 



(15) 



+ 



> CpCxCyS^ ^13 ■ (^^[S-p)^xSy + S^CxCySz 

S 13 

-T^ (c(5_o-)CxSy + C^Cj^Si'S^) — C23 • "So- (•S(5_o-)Ca;S2/ — S^CyS^Sz^ 



CxCyCz 
CyC^Sx / 



+C^Kyl |- (^^^ ~ ■ SpS^^j C1C2CXS2 - (^-^ - C23 ■ S^S^^j CiclsiSxj 

+Ciyl{(clsl - 4 {4 - slsl)) CzSz - (c^p+^^^CxS2 + C(^+^^)C2SiS^) C1C2 {4 " 4) 

+'2c^p-a+4>i-<P2f2CxCzSiS2SxSz] (16) 



I. .2 J (P-^) 



C12 

-C12 • C(^p-a) 



( 2 2 2\ ( 2 2\ 

%-<y) VV ~ - {(^(S+p-ct)(^x - C(S-p+a)^xJ 



2\ Cy^y^z 



(2 2 2\ ( 2 I 2\ 



CiiSySz 



Cx S:r 



, A-1 / \ ^p^y (2 2 2\ / 2 2\ 

+S13 ■ \ C{S-p)SxSy - CpCxCySzJ — — [C^ - S^S^J - S^S-p)~~ [^y ~ ^y) 



C-x^z 

^p^y 

^x^z 



+C13 • {s^5-p)SxSy + SpCxCySz) ^ (4 - S^S^) + C(5_p)^ - sfj 



/•-I / , A '^a'^y (2 2 2\ , '^x ( 2 2\ 

~S,23 ■ \ C{5-a)CxSy + C^CyS^Sz) — — \ - C^SzJ + S(5_o) — [Cy - SyJ 



+ 



C23 ■ {^(S-a)'^xSy — S^CySxS^ 
^^{<l>l-<l>2)^{p-a) 



^rj^y 



(2 2 2\ Cx / 2 2\ 

\Sx - C-x^z) - '^{S-a) — y^y - ^y) 



12 



+ C12 • C(p_^)S(^j_^2) 



C2S1S2 



^01 / Z' 2 _ 2 ^ 

A [^p'^y^y \^x ^x^n 
Us ^ ^ 



\S-p) 



{4 



Sy CxSxS. 



+C13 • {cpCySy (4 - 44) + C(S-p) {4 - 4) (^xSxSz) 



C1C2S2 



+ 



C23 



{s^CySy {sl. — 4^V) + *(<5-(t) — •Sy) CjjSjjSx) 
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CyC^SxSyS^ 



2 2\ ^x^^y^x^z 



+ 



(2(2 2\ . 2 2\ 



C1C2 



Cx^^z^x^ z 



+S(<,+^2)C2CxSl (clsl + sl (cl - sl)) 

/ 2 / 2 2\ 2\ /1 o 2 2M C1C2 

- S(5-p-,^i)"S2Sa: [C^ - \ Cy - Syj S^j -S (^s_^_^^^C2CxSi (^i - ZCyS 

_ 2 _ 2^ C2S1S2 



2c, 



C2CyS\S2Sz 



CyC^Sy 



(17) 



■ _ r^i 2 I '^(p-g) r / 2 2 2\ / 2 2\ 

P — '^kVt \ ~7 Y^{p-(j)(^xSx \ Sy - CyS^j - [C(^g_^_p_^^C^ - C(^g_p^^^S^j CySyS^ 



^(p-a)'^xSx [Sy — C,,S 



2 ,.2 „2 
y^z 



"5(5+P-(t)C^ + CySySz 



+ 



C12 
C12 ■ %-<j) 

-y- (c(S-p)^xSy — CpCxCyS^ + Cl3 ■ Cp {s(S-p)SxSy + SpC^CyS^ 
S,13 

S 

C23 



Sx 

Ccr. 



'^y i^x^z ^z) 



— {c(^S-a)^xSy + C^CySxS^ + C23 ■ C^ (s(^g_^^CxSy — S^CyS^S^ 



■^y \^x^z 

Cx^z 

2 ^ 

CyC^Sx I 



H<t>i-4'2f(p-'y) 



12 



+ C12 • C(p_a)S (</,!-</,; 



'2) 



C2SlS2S^ 



+ I ^T-^ + Cl3 • CpS^^ h — + C23 • C^S^^ 



C: 



13 



Cx.Cz^z 



,23 



C1C2S2 (1 - clsl) Cicjc^SiSx 



+ s 



{<T+4>2) ' 



C2S1S2SX 



{p-a+4>i-<l>2)' 



(18) 



■ _ /^i 2 I (p-g) r / 2 2 2\ / 2 2\ 

^ — '^kUt i ~ [(^{p-a)^xSx [Sy — CyS^J — (^C(5^^_^)Cj. — C(^g_p_^_^>^S^j CySyS^ 



C: 



12 



+ 



.Cl3 

S 

C23 



('^(5-p)'^a;'^J/ CpCxCySzj + Cl3 ' Cp {Sfg_p\SxSy + SpCxCySz 

L3 

— (c(5_g.')C-j;Sj^ + C^CySxS^ + C23C(j ^S(5_(^)C3;Sy — S^CySxSz^ 

'^(<?ii-<?i2)=^(p-a) , > ^ ^ \ C2C^SlS2 

h U2 • C(p_^)S(^^_^2) - 



2 



(C^s2 - sl) 



C: 



+ 1 + • c s^, + + C23 • c^S02 TT^ 



8 



I, '^Z^X^Z ^x ) 

An outline of the derivation is given in Appendix B. Concerning these equations, two 
remarks are in order: 



• As mentioned above, we can obtain exactly the same formulae as in Ref. [23] (and also 
in Ref. [27] but with a slightly different phase convention) if we set yi, — 0. This serves 
as a check of our derivation, at least for the part below the seesaw threshold. 

• p and a are only determined up to mr {n = 0, 1, 2, 3, ■ ■ ■) in Eqs.(6) and (9). Such an 
ambiguity in p and a also leads to an ambiguity in (pi and 02 defined in Eqs.(lO) and 
(11). However, both ambiguities cancel on the right hand side of Eqs.(14)-(19). So 
this ambiguity is harmless, as it should be. 

In addition to equations given above, a knowledge of the RGE evolution behavior of 6*1, 6*2, 4>i 
and 02 will be helpful in our following discussions. So for completeness, we have also 
derived one-loop RGEs for parameters in Ui, under the condition that eigenvalues of Y^, 
are hierarchical. We find that RG corrections to 61,62, 63, Si,, (f)i and 02 can be strongly 
enhanced by the factors Cij^ defined in Eq.(12), even if M3 and Aqut are only one or two 
orders apart in magnitude. The full analytical formulae are given in Appendix C. 



III. RG CORRECTION TO NEUTRINO PARAMETERS WITH THREE 
TYPICAL MASS PATTERNS: FROM Mz TO Agut 

We have verified our analytical formulae in Eqs.(13)-(19) by comparing their numerical 
solution with those obtained in the more conventional way, which is to integrate RGEs for 
Yi and k numerically, and then to calculate left-handed neutrino masses, leptonic mixing 
angles and CP- violating phases by diagonalizing these two matrices at different energy scales. 

In this section, we compare main features of the numerical result with those predicted 
in Eqs.(13)-(19). We try to clarify what corrections are possible for neutrino parameters, 
during the RG evolution from Mz to Aqut- Such a study is important in that, it helps us 
understand what values are possible or even preferred for neutrino parameters at the scale 
of Grand Unified theories. In the numerical calculation, we follow a bottom-up procedure. 
We start with the best fit values of neutrino parameters at the low energy, and numerically 
integrate Eqs.(13)-(19) to obtain left-handed neutrino masses, leptonic mixing angles and 
CP-violating phases at different high energy scales. Though such an approach may not seem 
well motivated from the perspective of a fundamental theory, it is quite advantageous for 
our present task. In this way, we no longer have to tune parameters at the high energy scale 
to meet low energy constraints. 

At the present time, low energy experiments have measured leptonic mixing angles and 
neutrino mass squared differences to a reasonable degree of accuracy (best fit values and 3a 
errors [33]): 
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9x ~ 33.2 ^4 go 5 9y ~ 45.0 i 9 3° j 13.0 5 

Am^i ^ml-ml^ (7.9 IJ;^) x lO'^eV^ ; 

^2.2 tli) X 10-^eV^ . (20) 



— mi 



However, we still lack a lot of information. Wc do not know about the smallest mixing angle 
9z and the absolute scale of neutrino masses, except for a few upper bounds. And we know 
nothing about leptonic CP-violating phases and the Yukawa coupling matrix Y^, at all. We 
can only speculate that Yi, might have some similarity to quark Yukawa coupling matrices. 
For such reasons, it is not practical to scan the whole parameter space for all possibilities. 

On the other hand, earlier works [17,18,27] have emphasized that enhancing factors like 
(^j^ play a significant role in the RG evolution of mixing angles and CP-violating phases, if 
left-handed neutrinos are nearly degenerate. So to be relevant and illustrative, we discuss in 
detail the RG correction in theories with three typically interesting neutrino mass patterns: 

(i) Normal Hierarchy: mi << m2 << ma. 
With Eq.(20), we find 



m2 ~ VAmii ~ O.OOQeV; mg ~ \J\Aml^\ ~ 0.047eV. (21) 

Then mi ^ 10~^eV is small enough to make neutrino masses hierarchical: 

Cfs' ^ -1.25, Cfs' ~ -1-0, C23' ~ -1-5 : = 10"^; 

« -1.0, Cfa^ ^ -1.0, (23^ ^ -1.5 : mi < 10"^. (22) 

Furthermore, from the determinant of Eq.(4) 



M,^^-f?.; ,.ir,r,f/('^M^) , (23) 



a 



where m„ = mctx [mi, m2, ma], while m^, and niy denoting the other two lighter left- 
handed neutrino masses. When ~ 1, we have in the most interesting case y^, ~ 0{1), 

1742 

Ma GeV ~ 6.4 x lO^^GeV. (24) 

ma 

(ii) Near Degeneracy: mi ^ m2 ^ ma. 

Left-handed neutrinos are nearly degenerate if the absolute mass scale is much larger 
than values given in Eq.(21). /3 decay, /3/3oz/ decay, and cosmological and astrophysical 
observations have all set upper bounds on certain combinations of left-handed neutrino 
masses [34] . A rather stringent upper bound on nearly degenerate left-handed neutrino 
masses is [35] 

^1,2,3 < 0.34eV. (25) 



10 



If left-handed neutrino masses lie rightly beneath this bound, we find that 

Cr/ ~ -5600, Cr/ « -201, C23 ~ -208 . (26) 

Prom Eq.(23), when ^ ~ 1 and ~ 0(1), 

1742 

M3 GeV ~ 8.9 X lO^^GeV ~ lO^^GeV. (27) 

ma 

Note that the pattern ma ^ mi ^ m2 is also possible, we shall consider this case in a 
different work. 

(iii) Inverted Hierarchy: ma << mi ^ m2. 
In this case 

mi ~ ^|Amii| ~ 0.047eV, ms f« ^fi^^\^A^^ ~ 0.048eV. (28) 

Then ^ 5 x 10~^eV is small enough to make a hierarchy between itself and the 
other two masses: 



cr/^-112, 


Cl3 ^ 


" 1-2, C23 ~ 


1.2 


: mg = 5 X 10"W; 


C12 ~ -111, 


^13 


' 1-0, C23 ~ 


1.0 


: ma < lO'^eV. 


Furthermore, when ~ 1 


and 


^O(l), 







(29) 



1742 

Ma GeV ~ 6.4 x lO^^GeV. (30) 

m2 



Note that in cases (i) and (iii), we shall always refer to the first lines of Eqs.(22) and (29) 
in our numerical calculation. And we shall also need the value of Ma (mostly in the factor 
In Ma) in our discussions. Since varying (In M3) by one order of magnitude requires changing 
the magnitude of Ma by a factor of C(IO^), the values of Ma given in Eqs.(24), (27) and (30) 
are "precise" enough for an order of magnitude estimation. However, it should be stressed 
that although we always assume yl ~ there can be significant errors in the subsequent 

discussions if ^ >> C(IO^) or ^ « 0(10"^) in Eq.(23). But such errors are easily corrected 
by taking into account the precise value of q. 

In the following subsections, we firstly discuss the evolution behavior of left-handed neu- 
trino masses and of enhancing factors (^^ defined in Eq.(12), then we study RG corrections 
to leptonic mixing angles and CP- violating phases in theories with each of the three neutrino 
mass patterns listed above. 
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A. RG Evolution of Neutrino Masses and Mass Ratios (, 



-1 



From Eq.(13), it is easy to obtain 



rrii (A, 



GUT J 



rrii (Mz 



exp 



t'=ln Agut 



IGtT^ Jto=lnMz 



Gi (t) dt 



1,2,3. 



(31) 



Since dominates (for i = 1,2,3) in general, the RG evolution of left-handed neutrino 
masses is mainly governed by a common scaling [25,27] 



m (Agut) 
rrii (Mz) 



exp 



1 /-i-i^ 
IGtt^ 7*0=111 



/:'=lnAGUT 



{t) dt 



1,2,3. 



(32) 



There is a sudden change in the direction of the common scaling at the point where i/i, ~ 0{1) 
is turned on, since y,y contributes to (t). This feature is quite obvious in Figure 1, where 
mi in cases (i) and (ii) and in case (iii) are plotted as functions of the energy scale, 
normalized by their values at Mz . 

Appreciable deviations from the common scaling may occur when Ut or y„ is large. Below 
the seesaw threshold, significant deviations are possible only in the MSSM when tan/3 is 
large, in which case i/r ~ C'(l). Above the seesaw threshold, appreciable deviations are 
generally possible since it is natural to have 7/,^ ~ (9(1). In the following, we shall elaborate 
on this problem in a new but more meaningful way, i.e. by studying the RG evolution of 
the factors (^'^ (for i < j, i,j = 1, 2, 3) defined in Eq.(12). 



1. Magnitudes of RG Corrections to 



As already mentioned, enhancing factors as Qj^ can be very important in the RG evo- 



lution of mixing angles and CP-violating phases. So it is necessary to study their own RG 
evolution behavior. From Eqs.(12) and (13), we find 



ttj — Co 



(33) 



Obviously, Qj would be constants if all left-handed neutrino masses varied by an exactly 
common scaling, i.e. = aj (for i,j — 1, 2, 3), which is of course not realistic. To estimate 
magnitudes of possible corrections to Qj (or (^^) from Mz to M3, we find that^ 



{Cij Cij) 



— 



-5. 



At 



(Cij Cij) 



Ciyl\n{M^/M^) 



167r2 



: in the SM 



- (C^^ - Cij) • (1 + tan^ 13) ■ 10-^ : in the MSSM 



(34) 

(35) 
(36) 



^Throughout this work, we shall omit the difference between Mz and the SUSY breaking scale 
for simplicity. 
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where we have used Eq.(33) and have taken ijr ~ rrir (Mz) / (174GeV) ~ 0.01 in the SM, 
and Vr ^ (Mz) / (cos/? ■ 174GeV) ^ 0.01/ cos/? in the MSSM. 

With Eqs.(22), (26) and (29), we find that vary httle from Mz to M3 in the SM, even 
in case (ii). If Eq.(26) is used, the variation of ^12^ is most significant: |A^]^Y^{^^| ~ 5.6 x 
10~^ , while variations of Cis^ and (23^ are much smaller: | ACi^^ / (i^^ \ ~ | A^:^^ / (23 I ~ 2 x 10~^ . 

In the MSSM, variations of are amplified when tan P is large. However, tan^ P alone 
is not large enough to generate appreciable variations of (^j^: \A(^^j^ /(^^'^\ ~ 2.5 x 10~^ in 
the case when tan/? ~ 50 but ~ 0{1). So corrections to in case (i) are always 

negligible from Mz to M3. Large variations of are possible only when \C^'^\ are also large 
enough. In the hmit \Cr/\ ^ 00, we find from Eq.(33) 

• ^ a,-aj a, - a, At CjyllnjM^/M^) 



(Mz) ^ 0, (Mz) + AO, (l+tan^/?)-10 
Q^Ms) C.,(Mz) " 0,(Mz) 



(38) 



When tan/? ~ 10, only in case (ii) is modified: \C,i2 (^z) /C12 (^3) I ~ 5.6. In con- 
trast, significant variations of ai'e popular when tan/? ~ 50. Now there are roughly 
\Ci2 (Mz) /Ci2 (Ms) I ~ 140, ICiV (Mz) /Cfs' (M3) I ~ (Mz) /C23' (M3) | - 5 in case (n), 
and |Ci2^(Mz) /Ci2^(M3) | ~ 2.5 in case (in). It is spectacular that |Ci2^| is diminished to 
orders of magnitude smaller at M3 than at Mz in case (ii). Numerically, we find that the 
typical magnitude of Ci2^ at M3 is of (9(100), while |Ci3^| and |C23^| are of a few tens''. 

Since the common scaling of ki (for i = 1, 2, 3) from Mz to Agut is only by a factor of 
0{1) and that Am^i (x kl — k\ = — ^2 {ki + /c2)^, the strong reduction in |^{^^| means that 
Am|i is orders of magnitude larger at M3 than at Mz. This point apphes also to |Am|2| 
(or |Am|;^|) when \C,23\ (or \C,i3\) is strongly damped. 

From M3 to Aqut, Uv ~ ^iX) usually dominates corrections to C,^^ both in the SM and 
in the MSSM with tan /? ~ 10, while Ut is comparable to in the MSSM when tan /? ~ 50. 
In whichever case, the dominant contribution is of the order 0{y,y), so we can make an 
estimation as in Eq.(38) 



5' (^3) ^ , ^ |/gln(AcuT/M3) 3x10- . 



1 + ^ _ ^ lo'^GeY. (39) 



(Agut) IOtt^ • (M3) 0, (M3) 

Obviously, corrections to in case (i) are again negligible from M3 to Agut- But it is 
remarkable that substantial corrections to Cjj^ are now possible in the SM, both in cases (ii) 
and (ni). With Eqs.(26) and (29) and also considering the fact (M3) | ^ (Mz) |, we 
find |Cr2^ (M3) /Cr2' (Agut) | - 170, Kf/ (M3) /Cf/ (Agut) | ~ |C2"3' (M3) /C2"3' (Agut) | ~ 6 in 
case (ii), and |Ci2^ (M3) /Ci2^ (Aqut) | ~ 3 in case (in). In the MSSM when tan/3 ~ 10, the 



^An interesting exception to this simple estimation is that ^23^ may vary little from Mz to M3. 
Such a behavior occurs when 9x is swiftly diminished to stay with 9z at very near zero, while 6y 
is kept at about 45° during the most part of the energy range from Mz to M3. One can easily 
understand this point with Eqs.(37) and (44). 
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situation is quite similar. In this case only of case (ii) is moderately damped from Mz 
to M3, but it still is a powerful enhancing factor in Eq.(39). For (i^^ and ^23^ of case (ii) and 
for (^12 of case (iii), the situation is the same as in the SM. In contrast, in the MSSM when 
tan/3 ~ 50, in general are strongly damped from Mz to M3. As a result, only Ci2^ of 



case (ii) may vary significantly from M3 to Agut, by the ratio |Ci2''^ (Ms) /(12 (Aqut) | ~ 3. 
However, there can also be appreciable modification in (^^ or ^2^^ of case (ii) or in of 
case (iii) if any of them is not so strongly damped from Mz to M3. 



2. Signs of RG Corrections to 
Now we turn to discuss the signs of corrections to By Eq.(37), we find 



0, (Ms) ^ 0, (Mz) + ^ 
Cij (Agut) ~ Cij (M3) + ^ 
where Ati = ln(M3/Mz), Atg = ln(AGUT/M3) and 



a, At 



ttj At 



167r2 ' 
2 



167r2 ' 







02 _ 




2 




Oi 




as _ 




2 




02 




03 _ 




2 





2 _ 2W 2 



r^l 2 \ 2 ( 2 2 2\ 22,0 

-C^Vr [Cy [C, - C^S^) - S^Sy + 2CsC^CyS^SyS, 

2 \ 2 ( 2 2 2\ 2 2 9 

''^kVt Y^y \^z ~ ^x^z ) ~ ^x^y ~ ^CgC^CyS^SyS^ 



x-y 



- C:yl {clsl - . 

- C:yl {clcl - si) ; 

- C:yl (4 si) cl ■ 



(40) 
(41) 

(42) 
(43) 
(44) 



In these equations, one should let y,^ = at energies below the seesaw threshold (or M3 
equivalently, in this work). Left-handed neutrino masses will keep their order of sequence if 
[tti — aj)/2 are of the same signs as {Mz) (for i < j; i,j = 1, 2, 3, respectively). 

From Mz to Ms, signs of (a^ — aj)/2 are determined by that of the C^y^ term. As 
estimated above, significant variations of are possible only in the MSSM, where Cl — l. 
It has been shown in Ref. [27] that tends to evolve toward zero whenever the correction 
is large; 6y is always smaller than but not far away from about 45°; and 6^ is always smaller 
than about 15°. By this we can conclude that the C^y^ term in Eqs.(42)-(44) are usually 
negative, and so that IC12I of cases (ii) and (iii) is enlarged with increasing energy scales in 
general. A positive correction to C12 (which diminishes IC12I) is possible only when Sz is large 
and is negative in Eq.(42), just as observed in Ref. [27]. However, it is numerically more 
difficult to arrange a positive correction to (12 in the whole energy range from Mz to M3 
than a negative one. For the same reason, C13 and C23 of case (ii) usually receive corrections 
of their own signs during the RG evolution. 

Prom M3 to Agut, the contribution from y^, is important both in the SM and in the 
MSSM. In the SM, the y,, term dominates (oj — aj)/2. Since 9i and 62 are totally arbitrary 
parameters, we can choose (a^ — aj)/2 to be negative or positive at M3 as we like. In the 
MSSM, the y^- term is also important when tan f3 is large. But we still can change signs of 
{tti — ttj)/ 2 by adjusting 61 and 62- Furthermore, mixing angles and CP- violating phases 
often vary dramatically at energies above M3 when \C,ij^\ are large. So (a^ — %)/2 may also 
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change their signs during the evolution from M3 to Aqut- In general, both positive and 
negative (a,; — o,j)/2 are possible at energy scales above the seesaw threshold. However, as 
will be clear in the following, signs of Qj are not likely to be changed despite this fact. 

In Figure 2, we illustrate the typical evolution behavior of {^—Cif^ = Ki^^l in the SM 
and the MSSM. It is obvious that (— CwM vary in a way just as discussed above. 



3. RG Evolution of in Nearly Singular Situations 



As explained above, signs of (oj — aj)/2 can be either negative or positive from M3 to 
Aqut; depending mostly on 61 and ^2- Immediately, there comes the question of whether it 
is possible to change signs of at energies above M3. We find that the answer is likely to 
be negative. 

When corrections to Qj are positive, are diminished and \Cij^\ niay be enhanced to 
extremely large values. Then the situation is nearly singular. However, we find that \Cij^\ 
are always dramatically diminished after they have reached some extremely large values. In 
such cases, extremely high but very narrow peaks appear in the plot of iCij^l against the 
energy scale. There seems to be a protecting mechanism that keeps ICij^l from going to 
infinity (or equivalently, keeps from going to zero). 

One can understand this point most easily in the SM, in which is small and so that 
Eqs.(42)-(44) are dominated by the C^^yl term. We also need Eqs.(C.l) and (C.2) in the 
limit — > 00: 

/ ^2 ^2 ^2 „2 ^2 ^2 \ 

cisi + • • •; (45) 

(46) 

It is crucial that, with (^^ given in Eqs.(26) and (29), the dominant correction to 62 is always 
negative, while 9i can either be enlarged by the ^12^ term or be reduced by terms of and 
(23'. As a result, 62 is always driven toward 0° when {Cri'l or \Ci^\ is large. In case (ii), since 
all three factors \Cij^\ are large, 9i can be driven toward either 0° or 90°, depending on the 
competition among different terms. In case (iii), since only \Ci2\ is significant, 9i is always 
driven toward 90°. 

Corrections to in case (i) are always negligible, so we only need to consider cases (ii) 
and (iii). We shall firstly discuss case (ii) in detail. 

In Eqs.(42), (43) and (44), a nearly singular situation is most likely to occur when any 
of 9i ~ 0° or 90° or 92 ~ 90° is satisfied: if there is only ^1 ~ 0°, a peak of |Ci2^| may appear; 
if there is only 9i ~ 90°, peaks of |Ci3^| and |C23^| niay appear; and if there is only 6*2 ~ 90°, 
peaks of |Ci2^|, ICis^l a-nd ICis] are all possible. But since in case (ii) 

Cr/ ^ -r-^ , (47) 

^,12 -r C,23 

a peak of |Ci3^| should always coexist with peaks of ICul and |C23^|- Such a situation is in 
fact quite rare. 
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Now wc can discuss how signs of dj are protected. 

Firstly, if (12 is driven to very near zero and ICul develops a peak, terms led by 
dominate Eqs.(45) and (46). Then 62 is dramatically driven to near 0°, while 6*1 is swiftly 
enhanced toward 90°. The smaller IC12I is, the more efficient this mechanism can take effect. 
As a result, (ai — a2)/2 in Eq.(42) is quickly driven to be negative, leaving no chance for (12 
to reach zero or even to get its sign changed. 

Secondly, if (23 is driven to very near zero and |C23^| develops a peak for some reasons, 
the (^23 term dominates Eq.(45). Then 9i is dramatically driven toward 0°. The smaller 
IC23I is, the more efficient this mechanism can take effect. Since 62 ~ 90° is not preferred 
in Eq.(46), c| 7^ in general. As a result, (02 — 03) /2 in Eq.(43) is quickly driven to be 
negative, leaving no chance for (23 to reach zero or even become positive. 

Finally, if |Ci3^| develops a peak (so there are also peaks of |Ci2^| and |C23^|), then ({2^ and 
Ci3^ in Eq.(46) can dramatically drive 6*2 to 0°. Since terms led by (12 C13 Eq.(45) are 
suppressed by s|, the (23 term then dominates Eq.(45) and drives 9i to 0°. Then (oj — aj)/2 
in Eqs.(42), (43) and (44) all become negative quickly. As a result, none of C12, Cis and C23 
will vanish or even get its sign changed. 

So it is not likely that signs of Qj can be changed. However, we should stress that such 
a discussion is only a way to understand how signs of Qj can always be preserved, while the 
latter point is in fact not proved. 

In case (iii), only \Ci2\ is important. A nearly singular situation is possible if 9i ~ 0° or 
62 ~ 90°. At the peak of |Ci2^|, 62 is dramatically driven to 0° and 9i is driven to 90°. Then 
there will be no more peaks. 

In the MSSM when tan/9 ~ 50, yr is comparable to y,y and the contribution from yr is 
not negligible in general. But numerically, we find that the mechanism discussed above still 
works as long as yi, is large. 

Numerical examples illustrating these features are given in Figure 3. 



B. The Normal Hierarchy Case 

When left-handed neutrino masses are hierarchical, iCjj^l ^t^'^ of ^(1)- So there are no 
enhancing factors in the evolution of mixing angles in Eqs.(14)-(16), while the evolution of 
CP- violating phases in Eqs.(17)-(19) can still be enhanced by the factor s~^. In the hmit 
Sz — > 0: 



+c:yi 



C: 



13 



+ Cl3C„S 



p*(5-p) 



C: 



23 



{5-a) 



+ Cl3CpS^J C^S2 + ( + C23%S^\ C2S1 



Sz 

C1C2 



C1C2 



+ 



(48) 



where we have neglected all terms that are not enhanced by s~^. It is remarkable that 
dominant corrections to S, p and a are exactly the same. It is also interesting that the term 
proportional to C^y^ vanishes in the limit C13 ~ (^23 ~ —1, while those led by yl become 
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proportional to (Cf' — C^). Since = Ci —1 in the MSSM, the contribution from y,, also 
vanishes in the limit C13 ~ C23 ~ — 1 in this case. 



1. RG Corrections in the SM 
In the SM, ~ rrir (Mz) /(174GeV) 0.01, and the typical contribution is 

^ In ^ (rad) ~ 10"^ (rad) ~ (lO'^') ° . (49) 

Sine we are mostly interested in cases with a large y^,, we shall consider only y,, ~ 0{1) in 
this work. When y,^ ~ (9(1), the contribution is 

yl . Agut, / 0.014 (rad) ~ 1.1° : M3 ~ lO^^Qgy 



167r2 Ms ^ ^ \ 0.029 (rad) ~ 1.9° : M3 ~ lO^^GeV ' 

Apart from these, all other miscellaneous terms (besides Cij^) on the right-hand side of 
Eqs.(14)-(19) usually damp the values given above strongly. Through out this work, we 
shall assume that the net effect of all these terms is equivalent to a factor of C(O.l). Though 
such an assumption is mainly based on our numerical experience and is far from precise, it 
serves as a crude estimation and can help us understand the most important part of RG 
corrections. With this assumption, we estimate that only corrections of (10^^)° are possible 
for mixing angles running from Mz to M3, while corrections of (9(0.1°) are possible from 
M3 to Agut, if Vu ~ 0{1). 

For CP-violating phases, significant corrections are possible in the energy range from M3 
to Agut, since there is the enhancing factor ^ 5. 

Firstly, the RG correction to 6z is negligible in the whole range from Mz to Agut, when 
9^ is of O (1° ~ 10°) at Mz. In this case, s;^ - C (5 - 50). From Eq.(50) and the O (0.1) 
factor explained above, corrections to GP- violating phases are of O (0.5° ~ 5°). 

Secondly, when 9^ is of 0(0.1° ~ 1°) at Mz, the RG correction to it is not neghgible 
from M3 to Agut- However, the order of magnitude of is not changed by the correction. 
So in this case, corrections to CP-violating phases are of O (5° ~ 50°), given that y^, ~ (9(1). 

Finally, though a value much smaller than possible radiative corrections may not seem 
natural for 9z, we shall consider the case 9^ « 0.1° for completeness. Since the radiative 
correction can enlarge 9z to O (0.1° ~ 1°) at energies above M3, corrections to CP-violating 
phases in this energy range are roughly the same as in the second case. But if 9^ is not 
magnified above M3, the corrections can be extraordinarily large. At energies below M3, 
corrections from yj. to CP- violating phases are still negligible when 9^ is of (0.01° ~ 0.1°). 
But when 9z < 0.01°, i.e. > 0.5 x 10^, the contribution from y^- can be so strongly 
enhanced as to become appreciable or even important. 

In Figure 4, we illustrate each of the possibilities discussed above. It is obvious that the 
phases S, p and a in the figure always vary by approximately the same size, just as predicted 
in Eq.(48). 



17 



2. RG Corrections in the MSSM when tan /? ~ 10 

In the MSSM when tan/3 10, ^ (Mz) /(cos/3 • 174GeV) ^ 0.1, and the typical 
contribution from which is 

^ In (rad) ~ 0.1°. (51) 

The contribution form is the same as given in Eq.(50). Then corrections to mixing angles 
arc of C(O.Or) in the range from Afz to Afg, but are still of 0(0.1°) from M3 to Aqut- For 
CP-violating phases, RG corrections from Mi to M3 can already be appreciable when 6z is 
of O (0.1° ~ 1°). This is in vast contrast to the SM case, where to make the contribution 
from appreciable, 9^ should be about one hundred times smaller than what is given here. 

Now we have a novel possibility that there can be large corrections to CP-violating phases 
both at energies above and helow M^. When 6z ^ O (0.01°) at Mz, the RG correction does not 
change the order of magnitude of 6z in the energy range from Mz to M^, so ~ (9 (5 x 10^) 
can enhance corrections to CP- violating phases to be of O (50°). At energies above M3 , 
may acquire a correction of C(0.1°) in general, but it is still of C(O.Ol) at M^. As a result, 
CP-violating phases can vary dramatically above M3, and corrections of (9(50°) are readily 
possible for them. In contrast, in the MSSM when tan/3 ~ 50, 6^ usually is magnified to 
be well above 0.1° in the energy range from Mz to M3. So corrections (from terms led by 
yl) to CP-violating phases are less enhanced and the phases vary httle at energies above 
M3. To conclude, it is only in the MSSM when tan/3 ~ 10 that large corrections to CP- 
violating phases are most probable, both at energies below and above M3. This observation 
is supported in Figure 5. 



3. RG Corrections in the MSSM when tan /? ~ 50 

In the MSSM when tan /3 ~ 50, yr ^ rrir {Mz) /(cos/3 - 174GeV) ^ 0.5, and the correction 
from which is 

Vr , AgUT 



The correction from y^y is still the same as given in Eq.(50). As a result, corrections of 0(0.3°) 
are possible for mixing angles running from Mz to M3. Prom M3 to Agut, corrections to 
mixing angles are still of C(0.1°). But for CP-violating phases, RE corrections can be 
extraordinarily large now. 

Firstly, the RG correction does not change the order of magnitude of when 6^ is of 
O (0.1° ~ 10°) at Mz. Then corrections to CP-violating phases are of C»(1.5°), 0(15°) and 
C(150°) when 9z ~ 0(10°), 0(1°) and 0(0.1°) respectively. When possible corrections are 
as large as 0(150°), CP- violating phases are often driven to near their (pseudo-) fixed points: 
the phases keep varying dramatically until the right hand side of Eq.(48) vanishes, and then 
they become stable against the energy scale. RG corrections to CP-violating phases can be 
largely damped because of this behavior. 
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Secondly, when 9z « 0.1° at Mz, the RG correction can enlarge 6z to be of O (0.1° ~ 1°) 
along the way from Afz to M3. However, extraordinarily large corrections to CP-violating 
phases can arise in a very narrow energy range at near Mz, if 9^ is extremely small at 
the beginning. In this case, (pseudo-) fixed points of CP-violating phase are often swiftly 
reached, and then the phases evolve little with the energy scale. Above M3, the phases 
start to vary again when the contribution from yy is turned on. However, the dominant 
contribution still comes from Ur- This is because 9^ is generally enlarged to O (0.1° ~ 1°) 
below M3, and so that the contribution from is less enhanced. 

Numerical examples illustrating these points are given in Figure 5. 

Apart from the magnitudes, signs of RG corrections to CP-violating phases are also 
important. In Eq.(48), the signs are determined by those of C^, C^, Cf, C13 and C23, and by 
values of CP-violating phases. It may not worth while discussing all possibilities, but one 
case is simple and interesting. In the SM, since C]^ = 1/2, Cj^ = —3/2 and since (±3 and (23 
are negative, terms led by y^ in Eq. (48) are always negative when both p + 4>i and a + (f)2 
are in the range (0 ~ tt), but the contribution from these terms is positive when the range 
is (tt ~ 27r). One can check this point with Figure 4. 

C. The Near Degeneracy Case 

When neutrino masses are nearly degenerate, corrections to mixing angles and CP- 
violating phases can be resonantly enhanced since enhancing factors \Cij'^\ in Eq.(26) are 
large. Before the discussion of any specific models, we remark that 

• In Eqs.(14)-(19), corrections to 9^, 6, p and a are enhanced by ^f^^, (^i^ and (^23 ^ while 
corrections to 9y and 9^ are only enhanced by and (^23 ■ 

• In Eqs.(14)-(19), the factors (^^ appear in terms led by yi, via three common com- 
binations: (12 ■ C2SiS2C(^j_<^2)) Cis"*^ ■ CiC2S2C^i and ^23"'^ ■ CiC^SiC^^. So there will be no 
resonantly enhanced contribution from y^, if 6*1 = 6*2 = (or 6*1 = 6*2 = 90° etc). 

• Just like c^^, C(^2 and C((^j_(^2) in the contribution from y^, certain functions of S,p 
and (7 can also be factorized out together with the enhancing factors Cjj^j both in 
contributions from i/^- and from y^. The result is collected in Table 1. With this, it 
is easy to find out what CP- violating phases can best damp the resonantly enhanced 
correction to a specific mixing angle or CP-violating phase. 

• In the contribution from y^, in Eqs.(14)-(19), the association of (^^ with CP- violating 
phases is rather simple. So it is easy to tell signs of corrections enhanced by different 
factors Cij^- We collect the result in Table 2, where we have assumed that minus signs 
in Eqs.(14)-(19) belong to phase factors. 

In the following we discuss RG corrections to mixing angles and CP-violating phases 
both in the SM and in the MSSM. 
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1. RG Corrections in the SM 



From Ml to M3, RG corrections are determined by the contribution from With 
the help of Eqs.(26) and (49) and also considering the overall O (0.1) factor explained be- 
low Eq.(50), we estimate that the correction (enhanced by Ci2^) to 9^ is of C(0.5°), and 
corrections (enhanced by and C^2-i) 9y and 6z are of O (0.02°). 

From M3 to Aqut, the contribution from ~ C^(l) is dominant and the magnitude 
is given in Eq.(50). Since C,^^{M^) pa C,^^{Mi), the correction (enhanced by (,^2) to Ox is 
of e> (1100°), and corrections (enhanced by Cii and (23^) to Oy and 9^ are of C (40°). The 
O (40°) correction to 9z is a quite generous gift: a correction of O (1°) could be interesting 
when 9z is really small at Afz [21,27,30]. Wc shall discuss this problem in more detail in 
the MSSM when tan/? ~ 50 (within this subsection, i.e. also for the near degeneracy case). 
For 9x1 however, the O (1100°) correction is overestimated. As we have discussed in the first 
subsection and also is vivid in Figure 2, in a large part of the energy range from M3 to 
Agut, |Ci2^I is more than an order of magnitude smaller than it is at M3. So the correction 
enhanced by (^12 should be about an order smaller in magnitude than estimated above. 

In cases when there are nearly singular situations, i.e. when \Cij^\ develop peaks, the 
"protective mechanism" discussed in the first subsection becomes important again in damp- 
ing extraordinarily enhanced corrections to mixing angles and CP-violating phases. As 
mentioned in the beginning of this subsection, corrections from xji, that can be enhanced 
by ai'c always controlled by Si and i-c Ci2^ by S1S2, Cis"^ by S2 and C^2Z by Si. 9i 

is dramatically driven to near zero whenever |C23^| develops a peak, and 92 is dramatically 
driven to near zero whenever \Ci2\ oi' \Ciz\ develops a peak. As a result, enhancing effect of 
all peaks of \C,ij^\ are strongly damped. However, the contribution from can still be very 
large in such a case (numerically, we find that a value of 0(100°) is popular). 

Signs of corrections to mixing angles are determined by CP-violating phases 5, p, a, 0i 
and 02 J and the competition among different contributions led by (^12^, ^13^ and (23- For 
an interesting example, we discuss in detail the case when all of the five phases arc in the 
range (0,7r/2). As discussed above, only the contribution from yi, at energies above M3 is 
important. 

• For 9 J.: In Eq.(14), the contributions from y,^ led by (12^ and Cis^ are positive, while 
that led by (23 is negative. Since ICu] is generally larger than ICis] and |C23^|, the net 
correction to 9x is positive in general, at near M3. However, since ^2 usually is more 
swiftly driven to near zero than 9i, negative contributing terms led by (23 Eq.(14) 
often have a chance to catch up with the positive contributions led by (^2 and (i^. As 
a result, sometimes 9x may make a change of its evolution direction and evolve toward 
zero. 

• For 9y: In Eq.(15), the part of the correction from y^, led by Cis^ is negative but that 
led by (2?^ is positive. Since |Ci3^| and |C23^| are comparable to each other in general, 
the two parts of the correction usually cancel each other to a large extent. So the 
net correction to 9y is relatively small. Furthermore, since 6*2 is often more swiftly 
diminished than 9i, the part led by ^23^ often is dominant and the correction to 9y is 
positive. 
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• For 6z- In Eq.(16), the correction from to 6^ is positive. So 9z only increases with 
the energy scale in this case. This feature leads to a quite interesting possibility that 
6*2 is comparable to the other two angles at Aqut, but it is diminished when the energy 
scale is decreased. 

Numerical examples illustrating these features are given in Figure 6. In the figure, 6y is 
decreased above M-^. This is because 6 > tt/2 in that energy range. 

For CP-violating phases, from Mz to M3, the factors \Cij'^\ alone can enhance corrections 
from Ut- to be of O (0.5°). Large corrections to CP- violating phases are possible only when 
9z is small enough. To make this point clear, we find (similar to Eq.(48)) in the hmit — > 0: 
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^C'^yl {\p+4>^)CxS2 + S(^+^^)C2SiS^) ^ + • • • , (53) 

where we have omitted terms that are not enhanced by s~^ (including those led by Ci2^)- 
The contribution from can be enhanced to O (0.02°) by C^^.^ and C^^, so ~ (9 (0.1°) is 
small enough to help generating corrections of C(10°). However, since \Ciz \ and 1^23^1 are 
comparable to each other in general, corrections from y^- are strongly damped if p ~ cr. 

From M3 to Aqut, in general is not an enhancing factor since a correction of O (40°) 
is possible to 6z- With C^^^ in Eq.(26) being the only enhancing factors, the situation for 
CP-violating phases is much like that for mixing angles discussed above. So in this energy 
range, corrections of O (100°) are possible for CP-violating phases. However, there is still a 
notable difference between the evolution of CP-violating phases and that of mixing angles. 
At energies right above M3, dz is still very small. So the enhancing effect of combined 
with that of Cj^^ and (^23^ in Eq.(53) can drive CP- violating phases to vary almost abruptly. 
Corrections of this origin are distinguishable since they are enhanced by s~^ and thus are 
exactly the same for 5, p and a. Only after 9^ has grown large enough, will the evolution of 
CP- violating phases slow down and acquire a strength comparable to that of mixing angles. 

The evolution of CP-violating phases is also illustrated in Figure 6. In the figure, the 
difference between the RG evolution of mixing angles and that of CP-violating phases (as 
discussed above) is obvious. 



2. RG Corrections in the MSSM when tan /3 ~ 10 

In this case, the contribution of is given in Eq.(51). From Mz to M3, C12 in Eq.(26) 
can enhance the correction to 9^ to be of O (50°), while corrections to 9y and 9z (enhanced 
by Ci3^ and (^23^) are of O (2°). Here, the 0(50°) correction to 9^ is usually negative and 
will be damped when 9,j. is near zero. A positive correction to 9,^. is also possible. But such 
a correction is always suppressed by factors of 0{sz)- Considering the upper bound for 
(< 1/5), a positive correction to 9^ should be smaller than 0(10°). The smaller 9z is, 
the smaller such a correction would be. If we take 9z '^'2,°, the correction should be smaller 
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than O (2°). For 6y, the situation is similar to that of 6x- The dominant terms arc usually 
negative, while a positive correction is possible but is always suppressed by s^. For 6^, the 
correction can be either positive or negative, depending on values of CP-violating phases 
and on the competition among the terms led by ({^ and ^23^- 

From M3 to Agut, Vu dominates the corrections and the situation for mixing angles is 
quite similar to that in the SM. A notable difference is that, now one can start the evolution 
with a small 6^ at M3, since it has a good chance of being damped from Mz to M3. 

For CP- violating phases, ^^2^ alone can enhance the corrections to O (50°), in the energy 
range from Mz to M3. But this correction is strongly damped when p ^ a. Alternative 
large corrections can come from terms led by Cis^ and ^23^, in cases when 6z is very small. 
The corresponding RGBs are the same as those in Eq.(53). Since dominant corrections to 
5, p and a are exactly the same, the relation p a can be easily retained in this case. CP- 
violating phases are often driven to near their (pseudo-) fixed points when 0^ is extremely 
small at Mz. 

From M3 to Aqut, corrections to CP- violating phases are dominated by and the 
situation is quite similar to that in the SM. 

For a numerical illustration and also as a check of the third remark at the beginning of 
this subsection, we stress that Ox and the phases 5, p and a cannot acquire their largest 
possible corrections simultaneously. For example, one needs p ~ cr to make the C,i2 term in 
Eq.(14) least suppressed, so that the correction to 9^ can be largest. However, in Table 1 
and Eqs.(17)-(19), the condition p a damps corrections (enhanced by Ci2^) to CP-violating 
phases most strongly. Also, this p a condition damps the RG correction to 9^ in the range 
from Mz to M3. One can understand this point with the help of Eq.(16) and the fact that 
Ci3^ and ^23^ are usually comparable to each other. These results are all illustrated in Figure 
7. Note that, since we take p f^i cr in the calculation, the correction to 9z is largely damped 
and 9z is kept unchanged until M3. As a result, CP- violating phases vary almost abruptly 
at the point where iji, is turned on, just as in the case of the SM. 

3. RG Corrections in the MSSM when tan /3 ~ 50 

In this case, the contribution of Hr is given in Eq.(52). Along with the usual 0(0. 1) factor 
explained below Eq.(50), Ci2^ in Eq.(26) can enhance the correction to 9x to C(1700°), in 
the energy range from Mz to M3. The other two factors Cis^ and C23^ in Eq.(26) can enhance 
corrections to 9y and 9z to O (60°). However, these values are overestimated. We have shown 
in the first subsection that, usually are strongly reduced from Mz to M3. In Figure 2, 
C,i2 is more than an order smaller in magnitude than it is at Mz , during a large part of the 
energy range from Mz to M3. So in general, the correction to 9x should also be an order of 
magnitude smaller than estimated above. The variation of \Cvi\ and l^^^^j is more moderate 
and the strength of the reduction is more uniform in the whole energy range than those of 
|Ci2^|. For such reasons, we re-estimate that the correction to 9x is roughly of O (170°), while 
those to 9y and 9^ are roughly of O (30°). 

Furthermore, (pseudo-) fixed points are always possible for mixing angles and CP- 
violating phases whenever corrections are very large. When certain angles or phases are 
near their fixed points, RG corrections to them are strongly damped. As a result, real 
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corrections to mixing angles and CP-violating phases depend not only on values estimated 
above, but also on how (pseudo-) fixed points can be reached. This is also true for RG 
corrections above M3. 

Prom M3 to Aqut, contributions from y^, and y,- are both important. As estimated in 
the first subsection, since |Ci2^| is generally orders of magnitude smaller at M3 than at Mz, 
and since ICis^l and |(^^^| are about 5 times smaller at M3 than at Mz, the correction to 6x 
is roughly of O (100°), while those to Oy and 0^ are of O (10°). 

• Por Ox'. Prom Mz to M^, the dominant contribution in Eq.(14) (that enhanced by 
but not suppressed by s^) is always negative. But if 9^ is large, there can be a positive 
correction. The correction to Ox becomes largest when |c(p_o-)| ~ 1. Since p = cr = is 
stable against RG corrections, this condition is easy to retain. Purthermore, |c(p_cr)| ~ 

1 can also lead to a large correction to in the energy range from M3 to Aqut, if 
there is in addition |c(<^^_02)| ~ 1- In contrast, if we want the correction to to be 
small, we need both C{^p^„) ~ and CpC^^ ~ c^-c^a ~ 0- These requirements can only 
be partially satisfied by 5 = p = 0i = (up to tt) and cr = 02 = 7r/2 (up to tt), which 
are also stable against RG corrections. 

• Por Oy-. Prom Mz to M3, dominant contributions in Eq.(15) (those enhanced by 
Ci3^) C-ii I but not suppressed by Sz) are also negative. Though a positive correc- 
tion seems possible in Eq.(15) when 6z is large, a numerical example is hard to find. 

The correction to Oy becomes largest when |c(5_p)| ~ |c(5-o-)| ~ 1 below M3, and 
\c(5-p)C(f,^\ ~ |c(5_o-)C02| ~ 1 (these two terms should be in opposite signs) above M3. 
These conditions can be satisfied, e.g. by5 = p = o' = 0i = O and (62 = tt, which are 
stable against RG corrections. In contrast, the correction is smallest both at energies 
above and below M3 when C(5_p) ~ C(^5-a) ^ 0. This condition can be satisfied, e.g. by 
5 = and p — a — /2, which are stable from Mz to Agut if there is 0i = 02 = 7r/2 
in addition. 

• Por Oz'. The sign of the dominant correction to Oz depends on CP- violating phases and 

on the competition among terms led by ^{^^ and (,2^ in Eq.(16). Prom Mz to M3, the 
O (30°) correction to Oz is quite spectacular: it means that a too small value is no longer 
natural for Oz- However, there is still a notable exception: if Cp ~ Cg. ~ in Eq.(16), 
the correction to Oz is strongly damped both at energies above and below M3, and so 
Oz can be kept at a small value. Only in this case, can a tiny Oz be probable. However, 
for this to happen, p and a must also be stable against RG corrections. In Eqs.(17)- 
(19) and Eqs.(C.4)-(C.6), we find that if S(p_o-) ~ (^(S-p) ~ C(^5-ct) ~ c<^i ~ c<^2 ~ 0) 
corrections to 5, p, cr, (pi and 02 can all be strongly damped, and so that 5, p, cr, 0i and 

02 are all stable against RG corrections. Up to tt, these conditions mean 5 — Q and 
p — a — (f)i — (f)2 — TT /2. In contrast, the correction to Oz can be large both at energies 
above and below M3, if CpC(^s-p) ~ ^c^c^s-a) and CpC^^ ~ CaC^2 > and all of them being 
of 0{1) in magnitude. A simple but interesting phase configuration that can satisfy 
these conditions is 5 = p + ir/i = a — ir/i = 7t/2 and 0i = 02 — tt = (which, however, 
is not stable against RG corrections). 

Por CP- violating phases, possible corrections are roughly of 0(100°) both at energies 
above and below M3, if only the factors are taken into account. But if the smallness of 
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9z is retained during a small energy range, there can be extraordinarily large corrections to 
CP-violating phases. As in previous cases, such corrections often drive CP-violating phases 
to near their (pseudo-) fixed points dramatically However, CP-violating phases usually 
need to be kept at special values if one wants to damp all large corrections to 6*^, just as we 
have mentioned above. So in such special cases, there should not be any large corrections 
to CP-violating phases, though 9^ may be tiny in a wide energy range such as from Mz to 

^GUT- 

• For 5: Eq.(17) for the running of 5 is rather complicated. But much simplified ap- 
proximation can be obtained in the limit 9z 0, which is given in Eq.(53). A no- 
table feature is that, only terms led by Cis^ and (^2^ ^"^^ possibly enhanced by s~^. 
The contribution from these terms can be dominant when 9^ is small enough [9^ < 

\Cl3/Cu'\-\C23/C{2'\-0{r)). 

• For p and a: Eqs.(18) and (19) for the running of p and a are also quite compli- 
cated. But the dominant contribution from iji, is now simple. We can predict signs 
of corrections enhanced by different factors with the help of Table 2, where the 
association of enhancing factors with CP-violating phases is clearly shown. Fur- 
thermore, much simplified approximations of Eqs.(18) and (19) can be obtained in the 
limit 9z ^ and the results are also the same as given in Eq.(53). 

In Figure 8, we illustrate cases in which the correction to a specific mixing angle is mostly 
enhanced or damped, just as discussed above. But the corrections shown are not largest in 
general. There can be larger corrections when specially chosen CP-violating phases are used. 
However, this is not our main concern here. What we want to demonstrate is that a good 
prediction of RG corrections can often be made with the help of Eqs.(14)-(19). For CP- 
violating phases, the situation is quite complicated and few general conclusions regarding 
their RG evolution can be reached. In Figure 8, we only give an example to illustrate how 
the RG evolution of CP- violating phases may be affected by 9^. 



D. The Inverted Hierarchy Case 

What is special with the inverted hierarchy case is that only Kul is moderately large 
among the factors \Cij^\ (for i < j; i,j — 1,2,3) defined in Eq.(12). So more interesting 
corrections can be expected than in the normal hierarchy case. However, since only |Ci2^| 
is large, the situation will not be so complicated as the near degeneracy case. Similar to 
Eq.(48), we find in the limit |^{^^| >> 1 and — >■ 0: 
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cl + O [s-') + • • • , (55) 

where the O (sj^) terms are exactly the same for 5, p and a, and have been given in Eq.(48). 
It is notable that terms led by Ci2^ are not enhanced by s~^. As a result, there are no more 
extraordinary corrections in this case than in the normal hierarchy case with extremely small 
Oz- However, the inclusion of a moderately large (^^2 leads to two non-trivial consequences: 
(a) the correction to Ox can now be much larger than in the normal hierarchy case, and (b) 
there can be much larger corrections to CP-violating phases than in the normal hierarchy 
case, when is not an efficient enhancing factor. 

1. RG Correction in the SM 

In the normal hierarchy case, enhancing factors ~ ^iX) and RG corrections to 

mixing angles are negligible (except for 9z when it is extremely tiny): (9((10^^)°) in the 
energy range from Mz to M3 and C (0.1°) from M3 to Agut- In the present case when \C,i2\ 
is large, only the correction to Ox is enhanced: O (0.01°) below M3 and O (10°) above M3. 
Note that the correction to 9^ becomes largest when \c{p-a)C{ij>i-(f)2)\ ~ 1 in Eq.(54). 

Corrections to CP- violating phases are dominated by (^^2 when iCi^"*^! >> (or equiv- 
alently, 6*^ >> IC12I ~ 0.01 ~ 0.5°). In this case, the situation for CP-violating phases is 
the same as that for 9x ■ So the contribution from is negligible and large corrections of 
C(10°) are possible only in the energy range from M3 to Aqut- In Eq.(55), contributions 
to CP- violating phases become largest when \c(^^^_^^)S(^p_a-)\ ^ 1- 

When 9z « 0.5°, corrections to CP-violating phases are dominated by and the 
situation is similar to the normal hierarchy case. 

In Figure 9, we illustrate the typical evolution behavior of 9x,9z,S,p and a in the SM. 
The competition between contributions from and is obvious. 

2. RG Correction in the MSSM 

In the MSSM when tan/? ~ 10, the contribution from Ur is about 100(~ tan^/9) times 
larger than that in the SM but the contribution from is the same. So corrections to 9y 
and 9^ are of O (0.01°) in the energy range from Mz to M3 and of O (0.1°) from M3 to Agut, 
while the correction to 9x (enhanced by Ci2^) is of O (1°) below M3 and of O (10°) above 
M3. Corrections to CP-violating phases are the same as that to 9x when \C,\2 \ » (oi" 
equivalently, 9z » 0.5°). 

In the MSSM when tan (3 ~ 50, the correction from yT is about 2500 times larger than 
that in the SM but the contribution from y^, is the same. So corrections to 9y and 9z are of 
O (0.1° ~ 1°) both at energies above and below M3, while the correction to 9x is of O (30°) 
in the range from Mz to M3 but is still of O (10°) from M3 to Agut- Note that in Eq.(54), 
the y^ correction to 9^ is always negative in the MSSM, while the sign of the contribution 
from yi, depends on simple phase factors. Corrections to CP-violating phases are the same 
as that to 9x when is dominant. 
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When 6z « 0.5°, corrections to CP-violating phases are dominated by and the 
situation is similar to the normal hierarchy case. This is true no matter tan/? is small or 
large. 

In Figure 10, we illustrate the typical evolution behavior of O^, 0^,6, p and a in the MSSM. 
In the calculation, we take p ~ o" to make the correction to 9x significant, and |p — cr| ~ 90° 
to make corrections to 9z,S, p and a significant. 

IV. SUMMARY 

In this work, wc have derived one-loop rcnormalization group equations for left-handed 
neutrino masses, leptonic mixing angles and CP- violating phases, both in the SM and the 
MSSM extended with three right-handed neutrinos. At energies above the seesaw threshold, 
we show explicitly the contribution from the Yukawa coupling matrix that connects left- 
and right-handed neutrinos. For simplicity, we have assumed hierarchical eigenvalues of 
this matrix in our derivation, so our analytical results may not be applicable when the 
eigenvalues are not hierarchical. And since we have also simplified the task by decoupling 
all right-handed neutrinos at a common scale, the discussion may have to be modified when 
the RG evolution between right-handed neutrino thresholds is important [15,28]. 

Based on these equations, we study possible RG corrections related to three typically in- 
teresting neutrino mass patterns: normal hierarchy, near degeneracy and inverted hierarchy. 

We firstly study the RG evolution of the factors Qj (for i < j; i,j = 1,2,3) defined in 
Eq.(12). We find that \Ci~/\ can be significantly damped from Mz to Aqut, both in the near 
degeneracy case and in the inverted hierarchy case. It is also possible that \Cij^\ m^'Y develop 
extremely high and narrow peaks, so that the situation is nearly singular. However, signs 
of (ij are not likely to be changed, neither is the order of sequence of left-handed neutrino 
masses. 

In the normal hierarchy case, RG corrections from Mz to Aqut are always negligible for 
mixing angles, except for 9^ when it is extremely small. Appreciable or even significant RG 
corrections to CP- violating phases are possible only when 6^ < 0{1°). In the SM, dominant 
RG corrections to CP- violating phases generally arise in the energy range from M3 to Aqut- 
In the MSSM when tan /? is large, dominant RG corrections generally arise from Mz to M3. 
Only in the MSSM when tan/3 is about 10, can large corrections to CP-violating phases 
arise both at energies above and below M3. 

In the near degeneracy case, possible large corrections to mixing angles and CP-violating 
phases are plethora. Mixing angles and CP-violating phases are often driven to near their 
(pseudo-) fixed points, since corrections are usually very large. Interesting mixing angles at 
high energy scales are often possible. For example, it is natural to find a large 9z (comparable 
to 9a:, 9y) at Agut- 

In the inverted hierarchy case, only and are significant enhancing factors. So the 
situation is much like that in the normal hierarchy case. However, because of the large (12 , 
the correction to 9x can be large, and significant RG corrections to CP-violating phases are 
possible even when 9~^ is not an efficient enhancing factor. 

To conclude, since RG corrections play a significant role in relating the low- and high- 
energy physics, an analytical understanding of the RG evolution behavior of neutrino param- 
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eters is necessary and important. Following earlier works, we have extended this understand- 
ing beyond the seesaw threshold by deriving RGEs for left-handed neutrino masses, leptonic 
mixing angles and CP-violating phases running at energies above the heaviest right-handed 
neutrino mass, under a few reasonable simplifications. The significance of these equations are 
demonstrated by studying the RG correction related to three especially interesting neutrino 
mass patterns. We expect that our work will be very useful for building realistic neutrino 
mass models at high energy scales. 
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Note added 

When this work is being completed, we notice that another paper about the similar topic 
is released by S. Antusch et al [36]. However, our strategy and approach are apparently 
different from theirs. The two works are complementary to each other. 



APPENDIX 

APPENDIX A: ONE-LOOP RGES FOR LEPTONIC YUKAWA COUPLING 

MATRICES 

Let Yu, Yd, Yi and denote Yukawa coupling matrices of up-type quarks, down-type 

quarks, charged leptons and the one that connects left- and right-handed neutrinos, respec- 
tively. At energies below the seesaw threshold, one-loop RGEs for Y/ and k are {Hf = YfYj 
for f = u,d, I, u): 

Yi^(ai + Ni)Yi- Ni^C\Hi- (A.l) 

« = + N^K + kNI; = &^Hi. (A.2) 

For coefficients, we have in the MSSM 

ai^-i^-gl + ^gtj+Tri'iHa + Hi)- = 3; 

a, = -{^-gl + Q9tj+Ti{m^); C^ = l; (A.3) 

and in the SM 

«i = - i^-^gl + \9i) + Tr (3iJ„ + 3iJd + Hi) ; C\=^-- 

«K = - (3^2 - a) + IV (6//„ + QHa + 2Hi) ; = -| (A.4) 
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At energies above the heaviest right-handed neutrino mass, one-loop RGEs for Yi,Yi, and 
Mr are 

Yi^{ai + Ni)Yi; Ni ^ C\Hi + H,; (A.5) 

Y,^{a, + N,)Y,; N^^ClHi + ClH,; (A.6) 

Mr = NrMr + MrNI- Nr = Cfl(yjy,)^. (aj) 



In the MSSM, 



In the SM, 



«i = - (^^?? + ^gi) + Tr (3/7, + ; (A.8) 

a. = - (^^? + 3^i) + Tr (3i/„ + H,) ; (A.9) 
= 3; Cr = 1; = 1; C,^ = 3; C« = 2. (A.IO) 



= - ( 7^?i + -A ) + Tr CiHu + 3//<i + //^ + ; (A.ll) 



4"" 4 
9 o 9 



q; 



. = - (^^i' + -/^ + Tr (3i/„ + 3//d + + H;) ■ (A.12) 

= I CT = -| = -| c^.^ = ^; ^« = 1- (A-13) 

If we make use of k also at energies above the seesaw threshold, we can derive from Eqs.(4), 
(A.6) and (A.7) 

k ^ a^K + N^K + kN'^ ; (A.14) 

a, = 2«, ; N, = K-Nr = CIHi + C';:H, , (A.15) 

where Nr = CrH^, , with Cr being the same as that in Eq.(A.7). Also, is the same as 
that in Eq.(A.2). In the MSSM, Q = 1; and in the SM, Q = 1/2. 



APPENDIX B: DERIVATION OF RGES FOR INDIVIDUAL PARAMETERS 

In the same way as in Refs. [18,27], we need to calculate Tmns = t^liNS^MNS to find out 
RGEs for leptonic mixing angles and CP- violating phases. Prom Eq.(8), 

where 

T„ = t/tf/, = -Ti ; Ti = UjUi = -TI . (B.2) 
Por T;, we find from Eqs.(A.14), (6), (7) and (10) 
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k' + T^k' + k'TJ = 2a^fi;' + TV^k' + k'TVJ ; (B.3) 
7V« = V\N^U, = Clul^sH'iUum + C^U^Kul . (B.4) 

In Eq.(B.3), for diagonal elements {i = 1, 2, 3) 

A;^ = 2 (q;« + ReN^.ii) h ; ImT^.^i = ImTV^.^^ = ; (B.5) 
and for off-diagonal elements (i, j = 1, 2, 3) 



k- — k- 

fi'i ~r nj 



(B.6) 



For Ti, we find from Eq.(A.5) 

Hi = 2aiHi + NiHi + HiNi . (B.7) 

Then from Eqs.(7) and (10) 

HI + ml + hit! = 2aiHl + NiH'i + H[Ni ; (B.8) 
Ni = U}NiUi = C///; + Cr (C/mnsC/.) (C/mnsC/.)^ ■ (B.9) 

In Eq.(B.8), for diagonal elements {i = 1,2,3) 

Hi,, ^2{ai + Ni.u) Hi,, , (B.IO) 
and for off-diagonal elements (i, j = 1, 2, 3) 



Note that Ti„ (for i = 1,2,3) are arbitrary, since Ui is only determined up to a diagonal 
phase matrix on its right. 

Furthermore, in order to derive RGEs for parameters in Y,^, we find from Eq.(A.6) 

H, = 2a,H, + N,H, + H,N, . (B.12) 

Then from Eq.(lO) 

HI + [% + Ulm,) Hi + Hi (t, + Ulm,)^ = 2a,Hl + N^Hl + HlN, ; (B.13) 
= Cl (t/MNst/.)^ Hi (UmnsU.) + C^Hl ; = Ultl, . (B.14) 

In Eq.(B.13), for diagonal elements (i = 1, 2, 3) 



29 



Hl,^ = 2 (a, + TV,,,) HI,, ; (B.15) 
and for off-diagonal elements (i, j = 1, 2, 3) 

= - fr^f/' N^., - (UIT^U.) . (B.16) 

u.ii u.jj ■' 

Just like the diagonal elements of Ti, Ty^u (for i = 1, 2, 3) are also arbitrary since t/, is only 
determined up to a diagonal phase matrix on its right. 

To calculate 9x, dy, Oz, S, p and & from Tj^^g , an auxiliary diagonal phase matrix is 
required on the left hand side of J/mns defined in Eq.(9), i.e. we have to use a more 
general parametrization of the MNS matrix in Eq.(B.l): 



^MNS — P - V ; P 



(B.17) 



where V denoting the original C/mns defined in Eq.(9). Then from Eq.(B.4) 

= ClV^H'iV + C^U^KUl . (B.18) 

Together with Eqs.(B.5) and (B.6), this equation shows that is independent of the phase 
matrix P. Furthermore, from Eqs.(B.9) and (B.ll), the product [P*TiP) is also independent 
of P. So in an equivalence of Eq.(B.l) 

yyt ^ VT^V'' - P* (Ti + P*P) P , (B.19) 

the equations of off-diagonal elements are obviously independent of the matrix P. They 
are all together six linearly independent equations of ^a;, 6y, 9z, S, p and a, and can thus 
determine these six quantities unambiguously. For the diagonal elements, {P*P),, {i — 
1, 2, 3) are only determined up to arbitrary (i = 1, 2, 3) , but this is of no problem since 
a, f3 and 7 are not physical by definition. We may choose whatever value for {i = 1, 2, 3) 
as we like in the calculation, or may simply ignore the equations for the diagonal elements 
in Eq.(B.19). 

In contrast, we can see from Eqs.(B.14) and (B.16) that to extract 9i, 62, ^3, S^, 0i 
and 02 from T^,, an auxiliary phase matrix on the right hand side of Ui, (just as the phase 
matrix P on the left hand side of J/mns) is not necessary. We can use the Ui, defined in 
Eq.(ll) directly during the calculation. There are totally six linearly independent equa- 
tions of 9i, 02, O3, (pi and 02 in Eq.(B.16), so these quantities can also be determined 
unambiguously. 



APPENDIX C: FULL EXPRESSIONS OF 61, 62, 63,6,,, (pi AND 02 

dF 

(Note that for what ever F : F = IQir'^ — ; t — \np, with p being the energy scale.) 

CLL 
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TABLES 



TABLE 1. Functions of S, p and a that can be factorized out together with (for 
i < j'l = 1)2,3) defined in Eq.(12), both in contributions from y^- and yi, in Eqs.(14)-(19). 
Note that terms led by Ci3^ and ^23^ in Eq.(17) are compHcated, so that the last two rows in the 
table do not apply to 6. 
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TABLE 2. Association of CP- violating phases with enhancing factors (for 
i < j'l hj = 1)2,3) defined in Eq.(12), in the contribution from yi, in Eqs.(14)-(19). Note that 
terms led by Ci3^ and ^23''^ in Eq.(17) are complicated, so that the last two rows in the table do not 
apply to S. 
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FIGURES 
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FIG. 1. The RG evolution of a left-handed neutrino mass (mi of cases (i) and (ii), and of 
case (iii)) from Mz to Aqut, both in the SM and in the MSSM. At Mz, 9^ = 33.2°, 0y = 45°, 
= 5°,S = 45° and p = C7 = 5°. At M3, = 0.8, = ^2 = 45° and 63 = = (f)! = (j)2 = 5°. 




FIG. 2. The typical evolution behavior of enhancing factors i^—Cij j defined in Eq.(12), both 
in the SM and in the MSSM. At Mz, 6^ = 33.2°, Oy = 45°, 9^ = 5°, 6 = 45° and p = a = 5°. At 
M3, = 0.8, 9i = 92 = 45° and 9^ = 6^ = (j)^ = = 5°. 
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FIG. 3. Typical peaks of enhancing factors \^—Cij J defined in Eq.(I2) in nearly singular situ- 
ations. At Mz, = 33.2°, 9y = 45°, 9, = 5°, S = 55° and p = a = 5°. At M3, = 0.8, Oi = 
e-2 = 89° and ^3 = (5^ = </>! = 02 = 5°. 



SM m ,,,^^ =115GeV 

e^(M^)-0.05° 
e^(MJ_0.0005° 




Log,„(u/GeV) 




8 
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SM m ... =115GeV 
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1 : 6^(M^)-5° 
2: 6^(M^)_0.05° 
3: 6,(M )-0.0005° 



Log,„(u/GeV) 



FIG. 4. The typical evolution behavior of the mixing angle Oz and CP-violating phases ((5, p, a) 
in the SM, in the normal hierarchy case. At Mz, 6x = 33.2°, 6*^ = 45°, 5 = 55° and p = (t = 45°. 
At Ma, = 0.8, Oi = 62 = (pi = h = 45° and 63 = 6^ = 5°. 
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FIG. 5. The typical evolution behavior of CP-violating phases {S, p, a) in the MSSM, in 
normal hierarchy case. At Mz, 9^ = 33.2°, By = 45°, 5 = 55° and p = a = 45°. At M3, = 
e^ = e2 = 4>i = <t>2 = 45° and e^ = 6y = 5°. 
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FIG. 6. The typical evolution behavior of mixing angles ( 
{d,p,a) in the SM, in the near degeneracy case. At Mz, Ox = 
and a = 5°. At M3, = 0.8, 6'i = 6*2 = 45° and (^1 = ^2 = O3 



'x,0y,9z) and CP-violating phases 



33.2°, 0y = 45°, 5 = 55°, p = 80° 
= 6^ = 5°. 
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FIG. 7. The typical evolution behavior of mixing angles {0x,0y,9z) and CP- violating phases 
p,cr) in the MSSM when tan/3 is small, in the near degeneracy case. At Mz, Ox = 33.2°, 
= 45°, 5 = 55° and p = C7 = 45°. At M3, = 0.8, 01 = 02 = 45° and c^i = = 03 = = 5°. 
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FIG. 8. The typical evolution behavior of mixing angles {0x,0y,9z) and CP-violating phases 
(5, p, a) in the MSSM when tan /3 is large, in the near degeneracy case. Figures A, B and C illustrate 
cases in which the correction to a specific mixing angle is mostly enhanced or damped. Figure A 
is devoted to 6x, but we have also included the other two angles in the plot. Similarly, Figure 
B corresponds to 6y, and Figure C corresponds to 6z- Figure D serves to illustrate how 6z may 
affect the RG evolution of CP-violating phases. At Mz, (9:,,. = 33.2°, 9y = 45° and 9^ = 0.05° 
(used in Figures A, B and C); at M3, t/u = 0.8, 9i = 62 = 45° and 63 = = 5°. Initial values of 
CP-violating phases (6, p and a at Mz, and 4>i and (p2 at M3) are marked on the graphs. 
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FIG. 9. The typical evolution behavior of mixing angles Ox and 9^ and CP-violating phases 



{5,p,a) in the SM, in the inverted hierarchy case. At Mz, 
and a = 5°; at M3, = 0.8, 6*1 = 6*2 = 45° and 4,1 = 4,2 = 



33.2°, 9y = 45°, 6 = 55°, p = 85° 
6. = 5°. 
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FIG. 10. The typical evolution behavior of mixing angles 9x and 6z and CP-violating phases 
{5,p,a) in the MSSM, in the inverted hierarchy case. At Mz, 0^ = 33.2°, Oy = 45° and 5 = 55°. 
Also at Mz, p = (T = 5° in the plot of 9x, but p = 85° and cr = 5° in the plot of 6z, 6, p and a. At 
Mg, = 0.8, 61 = 62 = 45° and (j)i = (f>2 = 63 = S^, = 5°. 
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